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Abstract 

We show the complete integrability and the existence of optical soli- 
tons of higher order nonlinear Schrodinger equation by inverse scattering 
method for a wide range of values of coefficients. This is achieved first 
by invoking a novel connection between the integrability of a nonlinear 
evolution equation and the dimensions of a family of matrix Lax pairs. It 
is shown that Lax pairs of different dimensions lead to the same evolution 
equation only with the coefficients of the terms in different integer ratios. 
Optical solitons, thus obtained by inverse scattering method, have been 
found by solving an n dimensional eigenvalue problem. 



The propagation of optical solitons in an optical fibre draws a lot of attention 
in recent times because of its plausible applications in telecommunication. It is 
known that the propagation of optical pulses of very short duration in a fibre 
is described by the dynamics of a nonlinear evolution equation, namely higher 
order Schrodinger equation (HNLS) [1-5]: 

d z E = i[ ai d TT E + a 2 \E\ 2 E] + e[a 3 d T rrE + a 4 d T (\E\ 2 E) + a 5 d T (\E\ 2 )E] (1) 

where, E represents an envelop of electric fields, z is the direction of propagation 
of optical pulse and r is time. The terms on the r.h.s. of ([!]) describe respectively 
the effects of group velocity dispersion (GVD), self-phase modulation (SMP), 
third order dispersion (TOD), self-steepening (SS) and self- frequency shifting 
via stimulated Raman scattering (SRS). 

If e = 0, (Q) reduces to nonlinear Schrodinger equation (NLS) and in this case 
optical soliton can be produced due to a balance between GVD and SPM. But 
the optical soliton thus produced, cannot propagate over a long distance || with 
their shapes undistorted, due to the propagation loss, as they travel along the 
optical fibre. This propagation loss in the fibre, however, can be compensated 
by utilising Raman effect, which is achieved by transmitting a pump light wave 
simultaneously through the fibre 0. It is interesting to note that the Raman 
process which takes care of propagation loss in the fibre, already exists within 
the spectrum of a soliton in HNLS equation (Q) [1,2,8]. It is the last term, which 
is responsible for stimulated Raman scattering. No additional pump light wave 
is needed to compensate the propagation loss for the ultrashort optical pulses 
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through the fibre. Moreover, for ultrashort pulses effects of higher order terms 
(last three terms) in (Q) cannot be neglected. It is these reasons the study 
of the soliton solutions of HNLS equation becomes an important issue for the 
propagation of ultrashort pulses in the fibre. But, unfortunately the complete 
integrability of HNLS equation for arbitrary coefficients o is not well understood. 
It is our aim, in this paper, to study the Lax integrability and consequently to 
obtain soliton solutions of HNLS equation (Q) by inverse scattering method 
(1ST) for a wide range of values of the coefficients o. Integrability of HNLS 
equation (1) by 1ST has so far been shown only for two fixed ratios of the 
coefficients [9,10]. 

Now, in order to obtain soliton solutions of HNLS equation by 1ST, it is 
convenient to study a gauge equivalent nonlinear evolution equation, 

d t q(x,t) = e[d xxx q(x,t) + "/i\q(x, t)\ 2 d x q(x, t) + j 2 d x {\q(x, t)\ 2 )q(x, <)] (2) 

which is related to the HNLS equation ([!]) cxplicitely through the following 
gauge transformation 

E(z,t) — exp[~ipx + iep 3 t]q(x,t) (3a) 
x = t + aipz (3b) 
t = a 3 z (3c) 

with p = oi/(3eo3) = 02/(604), 71 = 04/0:3, and 72 = 05/03. The equation 
is known as complex modified KdV (cmKdV) equation. Thus there is an 
one to one relation between the solitons of cmKdV equation and those of HNLS 
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equation through the relation (gja). Integrability of (|2|) is shown by many authors 
by various methods [9-13] only for two cases, when the coeffients of the terms 
on the r.h.s. are in the ratios (i) 1:6:0 and (ii) 1:6:3. Recently N solitary wave 
solutions of HNLS equation has been obtained for one parameter family of 7 
parameters ||. But the Painleve analysis restricts the integrability of HNLS 
equation further to a subset of 7 parameters, when the terms in the r.h.s. of 
are in the ratio 1:6:3. In this paper we show the complete integrability of 
HNLS equation by 1ST for all values of the 7 parameters. This is achieved by 
exploiting a novel connection between the dimensionality of the Lax pair and 
the integrability of HNLS equation. It is interesting to observe that the Lax 
pairs of different dimensions lead to the same equation of motion (Q), but 71 
and 72 in different ratios. Subsequently, we generalise the 1ST method for the 
n-dimensional Lax pair and obtain soliton solutions for HNLS equation for all 
integer valued ratios of the coefficients. 

In order to show the Lax integrability of HNLS equation vis a vis cmKdV 
equation, let us start with an n dimensional Lax equations, 

9 X * = U(x,t,\)V (4a) 
d t t> = V(x,t,A)* (4b) 

where, *&(x, t) is an n dimensional auxilliary field and the Lax operators U(x, t) 
and ~V(x,t) are n x n matrices of the form 

U = —iXH + A (5a) 
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V = eA xxx + e(A x A — AA X ) + 2eA 3 (5b) 
-2ieAS(A 2 - A x ) + 4eA 2 A - 4ieA 3 S 

In (^]) X is a c-no. diagonal matrix and the matrix A consists of dynamical 
fields, q(x, t) and q*{x,t) only. The explicit form of S and A may be given as 

n-1 

S = en - e nn (6a) 

i=l 

n—1 n — 1 

A(x,t) = '^a l (x 1 t)e in -^2a*(x,t)e ni (6b) 

i=l i=l 

where, eij is an n x n matrix whose only (zj)th. element is unity, the rest 
elements being zero and 0^(0;, t) represent the dynamical fields. They may be 
chosen either as q{x,t) or as q*(x,t). It is interesting to note that different 
choices of cti lead to the same cmKdV equation (|^) only with the coefficients of 
the last two terms in different ratios. To obtain the equation of motion for the 
dynamical fields q and q* , we first notice that 

S 2 = l, £A + A£ = (7) 

The equation of motion will then follow from the compatibility condition of (^) , 
namely 

d t U(x, t) - d x V(x, t) + [U(a>, t), V(x, t)] = 
and by using the relation (^) as 

ft A = ed xxx A - 3e(A 2 (d x A) + (d x A)A 2 ) (8) 

The above relation gives a nonlinear evolution equation for the matrix A. It 
is now straightforward to obtain evolution equations for q(x,t) or q*(x, t), by 



choosing an explicit form of A ([3t]). We will see the coefficients 71 and 72 
depend not only on the dimensions of the matrix A, but also on specific choices 
of the dynamical fields a,i as q or q* . Let us consider some specific examples to 
clarify this point. 

If we consider a two dimensional Lax operator, S and A in (^|) would be of 
the form 

£ = en - e 2 2 
= (1 

and 

A = a x e 12 - a\e 2 \ 

( ax \ 
{-a\ J 

and in this case only one dynamical field, a\ exists. Therefore, in two dimen- 
sional case the choice of Lax operator is unique, a\ = q and substituting A in 
(||), it is found that the equation of motion reduces to the well known Hirota 
equation ||, 

d t q = ed xxx q + 6e\q\ 2 d x q (9) 

where, the coefficients are in the ratio 1:6:0. 

However, for three dimensional Lax operator S and A are of the form, 

S = en + e 2 2 - e 33 A = aie X3 + a 2 e 2 3 ~ a*e3i ~ «2 e 32 
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We have, therefore, two possible choices of the Lax operators since two indepen- 
dent dynamical fields, ot\ and oli are present in A. In one case, we may choose 
both ai and a-i as q and consequently the equation of motion (|sj) becomes 

d t q = ed xxx q + I2e\q\ 2 d x q 

which reduces to Hirota equation @ by appropriate scaling of the fields, q and 
q*. In another possible case, we may choose a\ as q, but ai as q* and then 
both the Lax operator and the equation of motion become indentical to Sasa 
Satsuma case : 

d t q = ed xxx q + 6e\q\ 2 d x q + 3ed x (\q\ 2 )q 

where the coefficients are in the ratio 1:6:3. These two specific cases are already 
studied from various points of view. Moreover, no other ratios of the coefficients 
7 emerge upto three dimensions from the family of Lax operators, considered 
here. We will see, however, some interesting consequences as we consider four 
dimensional Lax operators. It follows from (^) that for four dimensional Lax 
operators, X and A takes the form 

S = en + e 22 + e 33 - e 44 

A = «iei4 + «2e24 + «3e34 — «ie4i — ck 2 e42 — «3e43 

and eventually, three different choices of the Lax operators are manifest in this 
case. One possibility, of course, may be considered when all the cti, and 
c*3 are chosen as q. The evolution equation (^) then, once again, reduces to 
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the known one, namely Hirota equation (f|) by appropriate scaling of the fields. 
Interestingly, however, if we choose both a\ and ct2 as q, but as as q* , the 
evolution equation (JsJ) becomes 

d t q = ed xxx q + 12e\q\ 2 d x q + 3ed x (\q\ 2 )q (10) 

where the coefficients are in the ratio 1:12:3. On the other hand, if we consider 
another possibility, where a\ is chosen as q, but both a-i and 013 as q*, the 
evolution equation turns out to be 

d t q = ed xxx q + 6e\q\ 2 d x q + 6ed x (\q\ 2 )q, (11) 

the coefficients being in the ratio 1:6:6. The last two cases are definitely new 
ones, which will be shown to be integrable. It is now clear that as we go to 
higher and higher dicmnsions, more and more new ratios of the coefficients 
appear for which cmKdV equation will be integrable. If fact, in general for an n 
dimensional Lax operator, as is evident from (^), [n — 1) number of dynamical 
fields, ai exist. We may, therefore, identify I number of a's as q and the rest 
(n — I — 1) ones as q*. Consequently, the evolution equation (||) yields as 

dtq = (d xxx q + 6le\q\ 2 d x q + 3(n - I - l)ed x {\q\ 2 )q, (12) 

which is nothing but cmKdV equation, where the coefficients are in the ratio 
l:6Z:3(n-/ - 1), i.c 

7i-6/, 72 =3(n-i-l). (13) 
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cmKdV equation is, thus, Lax integrable for all possible integer values of co- 
efficients, 71 and 72, which can be demonstrated by appropriately scaling the 
dynamical fields, q and q* in each case. It is evident now, for a given n dimen- 
sional matrix, (n — 1) Lax pairs can be constructed and each Lax pair of a given 
dimension, leads to (n — 1) different ratios of the coefficients. This is a significant 
achievement over the previous works. The existence of Lax pair, although, is 
a storng evidence for the integrability of cmKdV equation, it remains to show 
that the family of Lax pairs (|^) admit soliton solutions through 1ST. 

In order to obtain soliton solutions through 1ST in our case, we have gen- 
eralised the 3x3 AKNS type eigenvalue problem [10,14] tonxn eigenvalue 
problem. Notice that asymptotically the auxiliary field, *&(x,t) obeys a simple 
relation, viz. 



for the whole family of the Lax pairs (|5j) and it depends only on the dimensions 
of the matrix Lax pair. Consequently the scattering data matrix, which, by 
definition, connects the Jost functions, ^ l '(x, X)\ x =-oo to the Jost functions, 
ip( l >(x, \)\x=oo, f° r i — 1)2, •••,7i, [10,15] evolves with time in some universal 
form. Finally, solving n coupled Gelfand Levitan Marchenko equations one 
soliton solution may be expressed in a simple form 



t) = exp{{-i\x - 4ieA 3 i}S], 



(14) 



q(x,t) = {j]/y\n— l))sechA(x,t)exp(iB(x,t)) 



(15) 
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with 

A(x, t) = r\x — e(rj 3 — 3£, 2 T])t — 7 — —ln(n — 1) 
B(x,t) = £x + e(£ 3 - 3^ 2 )t + S 

where 7 and d are determined by the initial conditions and we assume that 
the simple pole is situated in the upper half plane at Ai = £ + iff} for one 
soliton solution. A detailed calculation of soliton solutions by 1ST method for 
an n dimensional eigenvalue problem is considered in |l5f| for a more general 
case than HNLS equation and will be published elsewhere. It is straightforward 
to obtain optical soliton in terms of the original parameters by substituting the 
relations (||) in ( p"5| ) and it turns out to be 

E(z, t) = {rj/yT(n - l))sechA(z, T)exp(iB(z, r)) (16) 

with 

A(z,t) = i]T - [ea 3 (j] 3 - 3£_ 2 f]) - T^—ri\z - 7 - \ln(n - 1) 

3ea3 3e z a 3 2"e^ai i 

It is interesting to observe from ( [l6|) that the envelop of one soliton solution 
admits a simple sech type shape, which can be easily produced by a mode- 
locked laser. Moreover, the intensity, I s = rf/(n — 1) and the width, T s = i]^ 1 
for one soliton in ( [l6|) are related by the expression I S Y 2 S — 1/ [n — 1). In terms 
of 7 parameters, given in (|l3|), the expression I S T 2 becomes I S Y 2 = 6/{37i + 
2 (72 — 7i ) } ! which is in full agreement with the result of |^] . 

9 



To conclude, the exact integrability of the cmKdV equation is shown for the 
integer valued ratios of the coefficients by establishing an intrcguing relationship 
of the coefficients 7 with the dimensionality of the Lax pairs. The family of Lax 
pairs, so obtained, is shown to admit soliton solutions. One soliton solutions 
are found by solving annxn eigenvalue problem through 1ST. For one soliton, 
the relation between the intensity and the width turns out to be simple and 
depends only on the diemnsions of the Lax pairs. 
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